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I present results for the two-loop soft anomalous dimensions for associated production 
of a single top quark with a W boson or a charged Higgs boson. The calculation uses 
expressions for the massive cusp anomalous dimension, which are presented in different 
'x" . forms, and it allows soft-gluon resummation at next-to-next-to-leading-logarithm (NNLL) 

accuracy. From the NNLL resummed cross section I derive approximate NNLO cross 
^ sections for bg — t- tW~ and bg — t- tH~ at LHC energies of 7, 10, and 14 TeV. 

^'- 1 Introduction 

The Large Hadron Collider (LHC) will produce top quarks via top-antitop pair or single top 
quark processes with relatively large cross sections. Given the importance of the top quark 
[1] to electroweak and Higgs physics, and the observation of single top events at the Tevatron 

*/^ ■ [2l m H] , it is crucial to have a good theoretical understanding of top quark production cross 

sections. An interesting channel to study is associated production of a top quark with a W 
boson, bg — )■ tW~ , which is sensitive to new physics and allows a direct measurement of the 

O I Vtb CKM matrix element. This process is very small at the Tevatron but has the second 

Q ■ highest cross section among single top processes at the LHC. A related process is associated 

production of a top quark with a charged Higgs boson, bg — > tH~ . Charged Higgs bosons 
appear in the Minimal Supersymmetric Standar Model (MSSM) and other two-Higgs-doublet 

r> I models (2HDM). In the MSSM there are two Higgs doublets, one giving mass to the up- type 

c^ I fermions and the other to the down-type fermions. Among the extra Higgs particles in the 

MSSM are two charged Higgs bosons, H~^ and H~, and the associated production of a top 
quark with a charged Higgs is a process that the LHC has good potential to observe. Since a 
central mission of the LHC is to find the Higgs boson and another is to look for supersymmetry, 
the associated production of a charged Higgs with a top quark is an important channel to study. 
The next-to-leading order (NLO) corrections to bg — )■ tW~ were calculated in [5] and to 
bg —7- tH~ in [6l[7l[8]. These processes are very similar with respect to QCD corrections and they 
have the same color structure. Soft-gluon emission is an important contributor to higher-order 
corrections, particularly near partonic threshold. The soft-gluon corrections can be formally 
resummed to all orders in perturbation theory. The resummation follows from the factorization 
of the cross section into a hard-scattering function H and a soft function S that describes 
noncollinear soft-gluon emission in the process [9l [10]. The renormalization group evolution 
of the soft function is controlled by a process-dependent soft anomalous dimension Ts- The 
calculation of Ts is performed in the eikonal approximation, which describes the emission of 



soft gluons from partons in the hard scattering and leads to modified Feynman rules in diagram 
calculations. At next-to-leading-logarithm (NLL) accuracy these corrections were resummed 
for tW~ production at the Tevatron and at the LHC in [TOl [HI 112], while the corrections 
for tH^ production were presented in [131 E]. These results involved the calculations of the 
one-loop soft anomalous dimension for these processes. 

Recent developments in two-loop calculations [151 IISI IH] have now made possible the re- 
summation of next-to-next-to-leading-logarithm (NNLL) corrections for QCD processes. Here 
we begin by calculating the two-loop soft (cusp) anomalous dimension for two massive quarks, 
and then using these results in the limit when one quark is massive (top quark) and one is 
massless (bottom quark) we calculate the diagrams for associated single top quark production. 
Since there are three colored partons in the partonic processes bg — t- tW~ and bg — )■ tH^ there 
are many diagrams to consider but the end result for the two-loop soft anomalous dimension 
for these processes can be written in a simple formula. We then use those results to calculate 
approximate next-to-next-to-leading order (NNLO) cross sections for tW~ and tH~ production 
at the LHC. 

2 Two-loop soft (cusp) anomalous dimension for a heavy 
quark- ant i quark pair 

We begin by presenting the calculation of the two-loop cusp anomalous dimension, which is the 
soft anomalous dimension for e~^e~ — )■ ti [151 IISI . 

We expand the soft (cusp) anomalous dimension as Ts = {as/'rc)Ts + {as/'rc)'^Ts + ■ • •, 
The one- loop soft anomalous dimension, F^ , can be read off the coefficient of the ultraviolet 
(UV) poles of the one- loop diagrams in Fig. 1. 

In the eikonal approximation, as the gluon momentum goes to zero, the quark-gluon vertex 
reduces to gsTp v^ /v ■ k, with gg the strong coupling, v a dimensionless velocity vector, k the 
gluon momentum, and Tp the generators of SU(3) in the fundamental representation. For 
example the integral for the diagram in Fig. 1(a) is given by 
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where i labels the quark and j the antiquark. The quark and antiquark velocity vectors obey 
the relations Vi ■ Vj = (1 + /3^)/2 and vf = t;| = (1 — /3^)/2, where 13 = Jl — A-m? / s with 
m the heavy quark mass and s the center-of-mass energy squared. The eikonal diagrams are 
calculated in dimensional regularization with n = 4 — e using Feynman gauge in momentum 
space. 

We find the one-loop soft (cusp) anomalous dimension 

(1 + /32), (l-p' 
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where Cp = {N^ — l)/(2iVc) with N^ = 3 the number of colors. This result can also be written 
in terms of the cusp angle [H] 7 = ln[(l + /3)/(l - /3)], with coth7 = {I + /3'^)/{2(3), as 

V'i^ = CF{lcoth-f-l). (2.3) 
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Figure 1: One- loop cusp diagrams with heavy-quark eikonal hues. 






(a) 



(b) 



(c) 






(d) 



(e) 



(f) 



Figure 2: Two-loop cusp vertex diagrams with heavy-quark eikonal lines. 



We now continue with the two-loop diagrams. In Fig. 2 we show graphs with vertex 
corrections and in Fig. 3 graphs with heavy-quark self-energy corrections. The dark blobs in 
Figs. 2(c) and 3(b) denote quark, gluon, and ghost loops. The black blob in Fig. 3(a) denotes 
three kinds of corrections shown in Fig. 4. We do not show graphs with gluon loops involving 
four-gluon vertices and graphs involving three-gluon vertices with all three gluons attaching 
to a single eikonal line since such graphs have vanishing contributions. For each diagram we 
include the appropriate one-loop counterterms for the divergent subdiagrams. The calculations 
are challenging because of the presence of the heavy quark mass. Dimensionally-regularized 
integrals needed in the calculation are shown in Appendix A. Using the results in Appendix A, 
the UV poles of the integrals for each diagram are provided in Appendix B. 

Combining the kinematic results in Appendix B with color and symmetry factors, the 
contribution of the diagrams in Figs. 2 and 3 to the two-loop soft (cusp) anomalous dimension 
is 
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Figure 3: Two-loop cusp self-energy diagrams with heavy-quark eikonal lines. 
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Figure 4: Detail of the black blob of Fig. 3(a) 
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where Ik denotes the integral for diagram k, e.g. I2d is the integral for diagram 2(d). Also 
l2cq and Jsfeg denote the quark-loop contribution in Figs. 2(c) and 3(b), respectively, while l2cg 
and I^bg denote the gluon-loop plus ghost-loop contributions to the respective diagrams. I^ai 
denotes the sum of the graphs 3(ali) and 3 (alii) detailed in Fig. 4 while I3a2 is the integral for 
the last graph in Fig. 4. On the right-hand side of Eq. (12. 4 p in addition to the two-loop soft 
anomalous dimension, F^ , which appears in the coefficient of the 1/e pole, there also appear 
terms from the exponentiation of the one-loop result and the running of the coupling which 
account for all the double poles of the graphs. Here f3o = (11/3)Ca — 2nf/3, with Ca = Nc 
and Hf the number of light quark flavors. From Eq. (12. 4p we solve for the two-loop soft (cusp) 
anomalous dimension: 
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where 
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We have written the two-loop result F^ in Eq. (12.50 in the form of a term which is a multiple 



of the one-loop soft anomalous dimension F^ , 



Eq. (12. 2p . plus a set of additional terms which 
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Figure 5: Leading-order diagrams for hg — )• tW . 



have been denoted as Mg. Here C2 = 7r^/6 and ^3 = 1.2020569 ■ ■ ■. The well-known two-loop 
constant K [19] is given hy K = Cyi(67/18 — (,2) — 5nf/9. The color structure of Tg involves 
only the factors CpCa and Cpnf. Note that as /3 — )■ 1, Mg — )■ (1 — (s)/'^- 
The result in Eq. (12.51) can be written in terms of the cusp angle 7 as 
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and is in agreement, but in a simpler and more explicit form, with the result for the cusp 
anomalous dimension of Ref. 



3 Two-loop soft anomalous dimension and NNLL resum- 
mation for hg — )► tW~ and hg — )► t}i~ 

We now turn our attention to processes that involve a bottom quark, a gluon, and a top quark 
as the colored particles in the hard scattering, namely tW~ and tli~ production. The leading- 
order diagrams for hg — )• tW~ are shown in Fig. [5l if one replaces the W~ by an li~ the graphs 
describe hg — )• tH^ . We treat the bottom quark as massless |13]. In this section we calculate 
the two-loop soft anomalous dimension that will allow us to resum the soft-gluon contributions 
to NNLL accuracy. 

In Fig. [6] we show the one-loop eikonal diagrams for these processes. Calculating the 
integrals associated with these diagrams we find the one-loop soft anomalous dimension for 
hg -> tW--. 
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where s = {pb+VgYi t = {Pb—PtY, u = {Pg—ptY, and mj is the top quark mass. The expression 
for hg — >■ tH~ is identical. This result is slightly different from the result in Ref. p^ [T3] because 
the axial gauge was used in those papers, while the result in Eq. (13. ip is calculated in Feynman 
gauge. Of course these differences are compensated by other terms in the resummed formula 
and the final result for the cross section is independent of the choice of gauge. 







Figure 6: One-loop eikonal diagrams with bottom quark- gluon-top quark vertex. 









Figure 7: Two-loop eikonal diagrams involving the bottom quark and gluon eikonal lines. 









Figure 8: Two-loop eikonal diagrams involving the bottom quark and top quark eikonal lines. 









Figure 9: Two- loop eikonal diagrams involving the gluon and top quark eikonal lines. 



To find the two-loop soft anomalous dimension we calculate the diagrams shown in Figs. U\ 
[HI m plus diagrams involving the top-quark self energy as in Figs. 3 and 4. Since they are three 
colored partons, with one of them massive, we calculate diagrams that contribute to the cusp 
anomalous dimension for each pair of partons using the results in the previous section in the 
limit when one or both partons are massless. Note that diagrams that involve gluons attached 
to all three eikonal lines either vanish or do not contribute to the two-loop result [20], and 
hence we do not show them. Combining the kinematic results for the integrals from Appendix 
B with color and symmetry factors we have 
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where h-ime denotes the terms involving gluons attached to all three lines that do not contribute 
at two loops and [hd\bt, for example, stands for the 2(d)-type diagram in Fig. 8 involving the 
b and t quarks. 

We thus find the two-loop soft-anomalous dimension for bg — > tW~ 
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where '^ s tw- ^^ given in Eq. (13. ip . The result for bg — ?■ tH is the same. 



With the two-loop soft-anomalous dimension at hand we are now ready to resum the soft- 
gluon corrections at NNLL accuracy. For tW~ production the resummed partonic cross section 
in moment space (with A^ the moment variable) is given by [9[ HOl [TTl [21] 
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and similarly for tH~ production. 

The first exponent |22l [23] in the above expression resums soft and collinear corrections 
from the incoming partons 
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where i stands for the incoming bottom quark {i = q) or the incoming gluon [i = g). Here 

Nq = N[{mly—u)/m^] and A^ = N[{m'^Y—t)/m^] where mw is the ly-boson mass. The quantity 



Ai has a perturbative expansion, Ai = Z]n(«s/7r)"Aj . Here A^^'' = Cp and A^^^ = CfK/2, 
while 41) = Ca and Af^ = CaK/2. 

Also A = En(tt./7r)"A^"\ with A^) = £>« = 0, and [1 



^.).C.C.(-|l.Hc,,^C3).C.,v(^-|) (3.6) 

andDf = (C^/Cp)Df. _ 

In the third exponent 7j/j is the moment-space anomalous dimension of the MS parton 
density 0j/j and it controls the factorization scale, fip, dependence of the cross section. We 
have ■ji/i = —Ai InNi + 7j where Ai was defined above, Ni = NiC^'^ with 7^; the Euler constant, 
and the parton anomalous dimension 7, = Z]n('^s/^)"7j where 7^"*^^ = "iCp/A and 7^^^ = /3o/4. 

jjbg^tw jg ^Yie hard-scattering function while 5'^'^"^*^ is the soft function describing non- 
collinear soft gluon emission [9l |T0] . The evolution of the soft function is controlled by the soft 
anomalous dimension Ts^tw-- Here N' = N{s/mf) with N = Ne'"^ . 

For tH~ production the resummed formula is essentially the same. The only difference, 
apart from the obvious use of the appropriate hard-scattering function for this process, is the 
definition of Ng and Ng. In this case, Ng = A^[(m^_ — u)/mjj-] and A^^, = N[{mj^- — t)/mjj.] 
where rriH- is the charged Higgs mass. 

The resummed cross section, Eq. (13.41) . can be expanded in the strong coupling, as, and 
inverted to momentum space, thus providing fixed-order results for the soft-gluon corrections. 
The NLO expansion of the resummed cross section after inversion to momentum space is 

a('^=a''^^^{csV,is,) + C2V,is,)} , (3.7) 

where a^ is the Born term for the process and /i/j is the renormalization scale. We use the 
notation Vk{s4) = [In (s4/m^)/s4]+ in tW~ production and Vi.{s4) = [In (s4/m^_)/s4]+ in 
tH~ production for the plus distributions involving logarithms of a kinematical variable S4 that 
measures distance from threshold (54 = at threshold). For bg — )■ tW~ , S4 = s+t+u—m^—m^, 
while for hg — )■ tH~ , Si = s + t + u — mf — m^- . The coefficient of the leading term is 



C3 



2(A«+AW). (3.^ 



The coefficient of the next-to-leading term, C2, can be written as C2 = d^ + ^2, with d^ denoting 
the terms involving logarithms of the scale and T2 denoting the scale-independent terms. For 
hg -^ tW- 

<« = -(4" + 4")i"(^) (3-9) 

and 

T, = -2Af i„(!!i-p) -2<) ,„(!!!|f^) -(^(.) + 4")iu(^) +2r<V- (3.10) 

For hg — )■ tH~ replace both my/ and irit in the above two equations by rriH-. 
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As discussed in [TOl [13] the expansion can also determine the terms involving logarithms 
of the factorization and renormalization scales in the coefficient, ci, of the 6{s4) terms. If we 
denote these terms as Cj*, then for tW~ production 
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while for hg — )■ tH~ replace both iriw and irtt in the above equation by niH-- The full virtual 
terms are not derivable from resummation, which addresses soft-gluon contributions, but can 
be taken from the complete NLO calculation. 

The NNLO expansion of the resummed cross section for hg — )■ tiy~ after inversion to 
momentum space is 
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For fe^f — )■ tif ^ again replace both mw and rrii by mj;/- in the above equation. It is important to 
note that all NNLO soft-gluon corrections are derived from the NNLL resummed cross section, 
i.e. the coefficients of all powers of logarithms in S4 are given in Eq. (I3.12p . from 2)3(34) down 
to Vq{s4). In Ref. [lU] and [13113, where NLL accuracy was attained, only the coefficients of 
1)3(54) and T>2{s4) were fully determined. Thus, at NNLL accuracy the theoretical improvement 
over NLL is significant. As discussed in [TOlD^ additional 6{s4) terms involving the factorization 
and renormalization scales are also computed. 

4 NNLO approximate cross sections for tW~ and tH~ 
production at the LHC 

We now use the results of the previous section to calculate approximate NNLO cross sections 
for bg -^ tW- and bg -^ tH' at the LHC. 

We begin with tW~ production. As has been shown in (TUl [H] the NLO expansion of the 
resummed cross section approximates well the complete NLO result for both Tevatron and LHC 
energies. In fact when damping factors are used to limit the soft-gluon contributions far away 
from threshold, as was also used for tt production [25] and s-channel single-top production [T7] . 
then the approximation is excellent. This shows that soft-gluon corrections are dominant for 
this process. 
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NNLO approx (NNLL) tW' cross section (pb) 


mt (GeV) 


LHC 7 TeV 


LHC 10 TeV 


LHC 14 TeV 


170 


8.24 


20.3 


43.6 


171 


8.09 


20.0 


43.0 


172 


7.94 


19.7 


42.4 


173 


7.80 


19.4 


41.8 


174 


7.66 


19.1 


41.2 


175 


7.53 


18.7 


40.6 



Table 1: The bg — )■ tW production cross section in pb in pp collisions at the LHC with 
y/S = 7 TeV, 10 TeV, and 14 TeV, with fi = mt and using the MSTW2008 NNLO pdf [26]. 
The approximate NNLO results are shown at NNLL accuracy. 



In Table 1 we provide numerical values for the tW~ cross section at the LHC for energies 
of 7, 10, and 14 TeV and a range of top quark masses from 170 to 175 GeV. The NNLO 
approximate corrections increase the NLO cross section by ~ 8%. We note that the cross 
section for bg — )■ iW^ is identical. 

At 7 TeV with rrit = 173 GeV the approximate NNLO cross section from NNLL resummation 
is 

^NNLOapprox^^^ = 173 GcV, 7 TcV) = 7.8 ± 0.2+°:^ pb . (4.1) 

The first uncertainty is from scale variation between mj/2 and 2mt and the second is from the 
MSTW2008 NNLO pdf at 90% C.L. At 10 TeV, again with m^ = 173 GeV, the cross section is 
19.4 ± 0.5t}:? pb, and at 14 TeV we find 41.8 ± 1.0+^^ pb. 

In Fig. [TO]we plot the bg — )■ tW~ NNLO approximate cross section from NNLL resummation 
at the LHC versus top quark mass for energies of 7, 10, and 14 TeV. 

Next we consider the process bg — )■ tH~ . The ratio of the vacuum expectation values, V2, f i 
for the two Higgs doublets is tan/3 = f2/fi, and the value of the cross section depends on the 
choice of this undetermined parameter. However, the overall percentage enhancement of the 
cross section from the higher-order soft-gluon corrections is independent of the value of tan/3. 

In Fig. [11] we plot the bg — )■ tH~ NNLO approximate cross section from NNLL resummation 
at the LHC versus charged Higgs mass for energies of 7, 10, and 14 TeV, using a value of 
tan/3 = 30. The NNLO approximate corrections increase the NLO cross section by ~ 15% 
to ~ 20% for the range of charged Higgs masses shown. We note that the cross section for 
bg —7- iH^ is identical (assuming the underlying model is CP conserving). 



5 Conclusion 

The cross sections for associated production of a single top quark with a W boson or with a 
charged Higgs boson receive large contributions from soft gluon corrections. These contributions 
were resummed in this paper to NNLL accuracy, thus extending previous NLL results. Attaining 
this accuracy requires the calculation of two-loop soft anomalous dimensions from the UV poles 
of dimensionally regularized integrals of two-loop eikonal diagrams. First the two-loop cusp 
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Figure 10: The cross section for tW production at the LHC with y/S = 7 TeV, 10 TeV, and 
14 TeV, and MSTW2008 NNLO pdf. 
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Figure 11: The cross section for tH production at the LHC with y/S = 7 TeV, 10 TeV, and 
14 TeV, and MSTW2008 NNLO pdf. 
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anomalous dimension was calculated, which is an essential ingredient to all NNLL resummations 
for QCD processes, and then the result was used to calculate the two-loop soft anomalous 
dimensions for bg — )■ tW~ and bg — )■ tH~. From the NNLL resummed formula approximate 
NNLO cross sections were derived and numerical predictions made for tW~ and tH~ production 
at LHC energies. These approximate NNLO corrections enhance the NLO cross section for tW~ 
production by ~ 8% and for tH^ production by ~ 15% to ~ 20%. 
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Appendix A: Dimensionally regularized eikonal integrals 

We list results for several dimensionally regularized integrals needed in the calculation of the 
two-loop soft anomalous dimension. 

"^"^ -(-l)-^-f vr^-f 2-^+^f r (l +'-)2Fr(il+'-;l; pA (A.l) 
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where 2-^1 is the Gauss hypergeometric function. 
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where Fi is the Appell hypergeometric function. 
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Appendix B: UV poles of the integrals for eikonal one- 
loop and two-loop diagrams for the soft (cusp) anomalous 
dimension 

Here we present the UV poles of the integrals for the one-loop eikonal diagrams in Fig. 1 and 
the two- loop eikonal diagrams in Figs. 2 and 3. 

First we list the integrals for the one-loop diagrams 
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Then we list the integrals for the two-loop diagrams: 






(B.l) 
(B.2) 

(B.3) 



/.^(i±£^iUln3^^-^^ 



S/32 



1 + /9 



In 



' 1-/3 ^ 
1 + /3. 



Li (^l^Il\ + C ' 



+ Li3 






(B.4) 



14 



hcq = nf 



(1 + /3') 



'2cg 



6/3 

5 (l + /3^) 
24 /3 



ri 


5 1 


[? 


6eJ 


ri 


31 1 


i?" 


30 ej 



In 



In 



• l-/3 \ 
U + /3, 



(B.5) 
(B.6) 



'2d 



(l + /3^) 
4/3 



^■"! 



/3^ 



+ 



-"(1^)- 



(3j e 
4/3 i 



In 



1 + /3/ 2 ll + /3^ 






'2e 



'2d 



(B.7) 
(B.8) 



'2/ 



2C2 + In 



2/1-/3 
U + /3. 



2/3 U + /3, 



12/3 I1 + /3, 



(B.9) 



'3c 



In terms of the cusp angle 
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the previous results for the integrals can be written as 
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The above expressions simphfy when one of the quarks is massless. In that case coth 7 = 1 
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The integrals hsted before then take simpler forms: 
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